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ABSTRACT 

A  description  is  given  of  some  experiments  made  to  examine  the  possible 
trapping  of  surface  waves  by  a  submerged  cylindrical  sill  in  an  otherwise 
uniform  ocean.  Theories  by  Longuet-Higgins  (1967)  and  Renardy  (1981)  indicate 
that,  at  certain  dangerous  frequencies/  nearly  resonant  conditions  may  obtain, 
leading  to  unusually  large  wavefields  over  the  sill. 

These  'resonances'  have  been  sought  experimentally  by  examining  two 
functionals  of  the  wavefield  over  the  sill  to  see  if  they  should  become 
unusually  large  near  certain  of  the  'dangerous'  frequencies.  The  results  show 
no  manifestation  whatsoever  of  the  resonances.  Possible  reasons  for  the 
discrepancies  are  considered. 

It  has,  however,  been  observed  that,  at  certain  locations  over  the  sill, 
the  waves  were  amplified  by  a  factor  of  four  or  five  times  their  magnitude  in 
the  ' ocean ' ,  irrespective  of  the  frequency  in  a  certain  interval .  This  could 
be  an  important  consideration  for  the  siting  of  rigs  on  prominent  topographic 
features  of  the  ocean  bed. 
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EXPERIMENTS  ON  WAVE  TRAPPING  BY  A  SUBMERGED  CYLINDRICAL  ISLAND 

+  t  *  t 

B.  J.  S.  Barnard  ,  W.  G.  Pritchard  '  and  D.  G.  Provis 

1.  INTRODUCTION 

The  purpose  of  the  present  paper  is  to  describe  some  experiments 
concerning  the  possible  trapping  of  waves  by  certain  topographic 
features  of  an  ocean  bed.  The  importance,  of  such  trapping  was  first 

•  i  *  * 

appreciated  by  Ursell  (1951,1952)  who  showed  that  the  Stokes  edge  wave 
over  a  plane  beach  (e.g.  see  Lamb  1932,  p  447)  corresponds  to  just  one 
of  a  set  of  discrete  eigensolutions  of  the  associated  mathematical 
problem.  The  amplitudes  of  the  wavemode s  corresponding  to  these 
discrete  eigenvalues  have  an  exponentially  decaying  structure  at  large 

i 

distances  from  the  shore  and  the  energy  contained  in  them  is  finite. 

For  this  reason  it  is  common  to  consider  the  waves  as  being  trapped  in 
the  zone  near  the  shore.  But  the  very  fact  that  the  waves  are  localized 
near  the  shore  mews  that  they  cannot  be  excited  directly  by  the  forcing 
effects  of  a  wavemaker  sited  at  a  large  distance  from  the  shore.  Thus, 
to  realize  these  modes  in  the  laboratory  Ursell  (1952)  forced  the  motions 
in  the  zone  near  the  shore  by  oscillating  the  side  walls  of  his  channel. 
He  also  suggested  that  it  might  be  possible  to  generate  them  thro’igh  a 
nonlinear  mechanism,  and  this  was  subsequently  shown  to  be  the  case  by 
Galvin  (1965). 

For  the  mathematical  problem  considered  by  Ursell  (1951,  1952)  the 
eigenvalues  are  all  real  but,  when  slightly  more  complicated  geometric 
situations  are  examined,  the  operators  associated  with  the  mathematical 
problem  often  are  not  self  adjoint  (cf .  Lozano  a  Meyer  1976) .  Thus  we 
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can.  In  general,  expect  any  discrete  eigenvalues  arising  In  these 
situations  to  be  complex,  the  imaginary  part  corresponding  to  a  'leakage* 
or  radiation  of  energy  to  the  far  field.  Nevertheless,  it  might  be 
anticipated  that,  under  certain  conditions  (in  particular  in  situations 
bearing  similarities  to  those  in  which  edge  waves  are  possible),  the 
imaginary  parts  of  some  of  the  discrete  eigenvalues  would  be  small 
and  we  could  expect  such  waves  to  have  properties  closely  related  to 
those  of  the  perfectly  resonant  situation.  This  has  indeed  proved 
to  be  the  case  in  some  specific  problems.  For  example,  it  .has  been 
shown  for  a  submerged  seamount  (Longuat-ftiggins  1967,  Renardy  1981) 
and  for  conical  islands  with  a  shore  (Shan,  Meyer  a  Keller  1968, 

Meyer  1971,  Lozano  l  Meyer  1976)  that  there  are  discrete  eigenvalues 
with  very  small  leakages.  The  fact  that  these  wavemodes  can  potentially 
have  a  large  response  in  certain  regions  of  the  ocean  (near  the 
important  topographical  features)  and  yet  'leak*  energy  to  the  far  field 
is  of  considerable  interest.  For,  in  contrast  to  the  edge-wave 
situation,  the  leakage  allows  for  the  possibility  of  exciting  these 
'nearly  resonant*  nodes  directly  (without  the  need  for  a  nonl  inear 
coupling  mechanism)  by  means  of  a  wavetrain  generated  at  large 
from  the  island. 

Only  two  serious  attempts  appear  to  have  been  made  so  far  to 
investigate  experimentally  this  concept  of  wave  'resonance*  in  open 
bodies  of  water.  One  of  these,  a  forerunner  of  the  present  paper,  was 
a  study  by  Pro  vis  (1975)  (and  see  Provis  1976)  who,  for  the  cases  of 
a  conical  island  and  a  submerged  cylindrical  sill,  considered  the 
motions  arising  near  these  topographic  features  from  the  action  of  a 
wave  maker  at  one  side  of  a  square  tank.  The  other  study  (Pita  1977) 
concerned  the  motions  above  a  submerged  cylindrical  sill.  Even  though 
both  these  studies  were  undertaken  with  considerable  care,  neither 
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experiment  produced  any  prlma  facia  evidence  to  suggest  the  presence 
of  resonant-dike  resonse.  Pite  aeasured  wave  amplitudes  along  a  radial 
line  normal  to  the  direction  of  propagation  of  the  incident  wavetrain 
and  found*  for  three  different  frequencies,  that  his  results  vers 
described  fairly  well  by  a  theoretical  model  incorporating  an  allowance 
for  the  effects  of  daaping.  It  was  inferred  from  this  agreement  that, 
under  suitable  conditions,  the  near  resonances  predicted  by  Longuet- 
Biggins  (1967)  should  be  evident.  But  Renardy  (1981)  has  shown  that 
the  frictional  effects  were  not  correctly  incorporated  into  the  theory 
used  by  Pite,  and  tbs  agreement  may  therefore  have  been  fortuitous. 

One  complication  of  experimental  work  in  this  area  is  the 
recognition  of  a  'resonance*  should  it  occur.  So,  in  order  to 

establish  a  procedural  basis  on  which  to  locate  'resonances',  we 

» 

decided  to  measure  mama  functionals  associated  with  the  wave  field, 

(certain  real -valued  properties  of  the  wavefield  as  a  function  of 

frequency, )  and  then  to  count  as  resonance  any  range  of  frequencies 

0 

for  which  a  functional  assumes  abnormally  large  values.  For  example, 
if  we  were  to  examine  the  maxima  wave  amplitude  over  the  entire  ocean 
(relative  to  that  of  the  incident  wavetrain) ,  we  should  refer  to  any 
frequency  for  which  this  quantity  takes  on  unusually  large  values  as 
being  a  'resonant'  frequency.  Not  only  are  these  functionals  useful 
in  characterising  any  resonances,  they  also  provide  a  convenient  way  of 
summarizing  the  results  for  comparison  with  theory,  as  we  shall  indicate 
below. 

The  particular  experiment  to  be  described  in  this  paper  concerns 
the  possible  trapping  of  waves  by  a  submerged  cylindrical  sill  in  an 
open  ocean  of  uniform  depth.  The  initial  work  of  Provis  (1976) ,  in  which 
trapping  near  a  conical  island  with  a  shore  was  considered,  did  not  reveal 
the  predicted  resonances.  A  possible  explanation  of  this  result  is 
that  the  resonances  were  overwhelmed  by  dissipation,  which  should  have 
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been  substantial f  in  ths  sons  near  the  shore  of  the  island,  (cf.  Mahony 
a  Pritchard  I960) .  So,  not  only  did  the  submerged  sill  offer  the 
attraction  of  being  relatively  insensitive  to  dissipative  effects 
(a  feature  that  has  subsequently  been  confirmed  in  a  calculation  by 
Renardy  1981)  but  also  the  explicit  solutions  given  by  the  theory 
of  Longuet-Higglas  (1967)  provided  a  framework  for  comprehending  the 
experimental  results.  He  therefore-  made  measurements  of  both  the 

wave  amplitude,  as  a  function  of  frequency,  and  of  the  'energy* 
functional.  There  was  no  evidence  from  these  results  to  suggest  the 
occurrence  of  the  dramatic  resonances  predicted  theoretically. 

Xn  an  attempt  to  find  a  possible  reason  for  this  disparity, 

» 

Renardy  (1981)  relaxed  the  shallow-water  assumption  used  by  Longuet- 
Higgins  and  established  a  complete  eigenfunction  expansion  for  the 
(linear)  problem.  These  more  exact  calculations  confirmed  the  'resonances' 
of  the  Longuat-Higgina  theory,  indicating  even  larger  responses  over  a 
slightly  narrower  bandwidth,’  but  gave  the  surprising  result  that  the 
locations  of  the  resonant  frequencies  were  shifted  substantially  from 
those  predicted  by  the  shallow-water  model.  Renardy  also  showed  that 
dissipative  effects  cm  the  plateau  section  of  the  sea  mount  should  have 
had  only  a  very  minor  influence  on  the  wave  amplitudes.  Thus  again  the 
empirical  results  were  at  variance  with  the  theoretical  predictions. 

The  reason  why  there  is  such  a  large  disparity  between  our 
experimental  results  and  the  theoretical  model  is  not  understood.  As 
is  the  case  for  all  laboratory  experiments,  there  are  several  respects  in 
which  the  experimental  conditions  did  not  match  properly  the  mathematical 
problem  studied  by  Renardy  (1981)  and  it  is  conceivable  that  any  of 
these  could  have  played  an  important  modifying  role  in  the  structure 
of  the  wavefield.  He  have  therefore  given  a  fairly  detailed  description 
of  the  experimental  apparatus  and  a  table  suamarising  the  main  experimental 
results  for  the  two  functionals  that  have  been  measured. 


The  tank  used  for  these  experiments  was  constructed  by  placing 
a  large  sheet  of  polythene  on  the  floor  of  the  laboratory  and  raising 
the  edges  of  the  sheet  over  concrete  walls.  A  plan  view  of  the  tank# 
together  with  its  dimensions,  is  given  in  figure  1.  Located  near 
three  sides  of  the  tank  were  plane  beaches  of  slope  approximately 
Isle.  These  beaches  each  sloped  towards  the  cylindrical  sill  near 
the  centre  of  the  tank.  A  large  proportion  of  the  fourth  side  of 
the  retaining  wall  of  the  tank  consisted  of  a  plank  of  laminated  tinkler 
This  plank#  which  was  used  to  generate  waves  in  the  tank#  was  hinged 
near  the  floor;  at  the  centre#  near  the  top  of  the  plank,  it  was 
attached  to  a  mechanism  used  to  force  the  wavemotion#  as  described 
below.  With  the  plank  placed  behind  the  polythene,  sheet#  no  leakage 
around  the  ends  of  the  wavemaker  was  possible. 


Located  near  the  centre  of  the  tank  was  an  upright  'cylinder'# 
referred  to  in  figure  1  as  a  sill#  of  mean  diameter  100.4  cm  and  height 
13.8  cm.  The  sill  consisted  of  a  disk  of  perspex  with  strengthening 
webbing  on  its  underside,  and  a  circumferential  skirt  extending  down 
to  approximately  5  mm  from  the  floor  of  the  tank.  The  structure  was 


mounted  on  three  feet  the  heights  of  which  were  adjustable  to  allow 
for  levelling  of  the  upper  surface.  The  levelling  was  done  with  the 
tank  filled  with  water  as  follows.  A  finely  pointed  cone  was  attached 
to  the  end  of  a  micrometer  screw  and  the  micrometer  was  clamped  to  a 
small  tripod.  The  tripod  was  then  placed  on  the  upper  surface  of  the 
sill  and  the  micrometer  was  slowly  screwed  down  until  the  tip  of  the 
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Figure  1.  (a)  A  plan  view  of  the  tank  used  for  the  experiments. 

(b)  A  sketch  of  the  sill.  (All  dimensions  are  in  metres.) 


cone  just  broke  the  water  surface,  the  micrometer  setting  at  this 
height  providing  a  measure  of  the  distance  between  the  upper  surface 
of  the  sill  and  the  water  surface.  After  the  sill  had  been  suitably 
levelled  a  systematic  survey  indicated  that  its  upper  surface  deviated 
from  a  horizontal  plane  by  no  more  than  0.2  mm.  This  'depth  gauge' 
was  also  used  to  set  the  water  level  before  each  experiment:  the 
tripod  was  placed  at  a  specific  location  on  the  sill  and  the  water  level 
was  adjusted  until  it  was  within  0.02  mm  of  a  nominal  value. 

Because  of  the  large  surface  area  involved  in  this  experiment, 
contaminants  could  easily  accumulate  at  the  water  surface  (originating 
both  in  the  air  and  the  water  itself)  and  unless  special  precautions 
were  taken,  the  resulting  surface  films  could  give  unacceptable  levels 
of  wave  damping.  He  did  mat  make  direct  quantitative  tests  of  this 
feature,  but  we  were  guided  by  our  experiences  from  other  surface-wave 
experiments  carried  out  in  the  same  laboratory  (see,  for  example,  Barnard 
ft  Pritchard,  1972,  Barnard,  Mahony  ft  Pritchard  1977,  Mahony  a  Pritchard  1980, 
and  Bona,  Pritchard  ft  Scott,  1981).  This  experience  suggests  that  a 
skimming  of  the  surface  just  before  an  experiment  leaves  it  clean  enough 
for  about  an  hour  before  the  build  up  of  contaminants  causes  significant 
changes  in  the  wave  motions.  The  skimming  of  the  surface  films  was  brought 
about  by  allowing  a  small  overflow  into  a  circular  weir  placed  in  the  tank; 
the  outflow  over  the  weir,  drawing  mainly  from  the  s  turf  ace  layers,  pulled 
off  the  surface  films.  This  clewing  procedure  ran  continually  whenever 
an  experiment  was  not  in  progress.  Then,  just  before  starting  a  new 
set  of  measurements,  the  weir  was  raised  and  the  flux  of  water  into  the  tank 
was  cut  off.  A  reference  measurement  of  the  wavefield  was  made  at  the 
beginning  and  end  of  each  experiment  to  establish  the  'drift'  of  the 
wavefield;  this  was  less  than  5%  for  all  the  results  given  here. 
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The  beaches 

ha  Indicated  above,  three  sides  of  the  tank  were  flanked  by 
beaches  extending  from  just  above  the  water  level  at  the  tank  wall 
and  sloping  at  lilO  down  to  the  floor.  The  beaches  were  made  of 
alusinium  and  were  covered  with  a  coat  of  priming  paint  which  was  rough 
enough  to  retain  a  film  of  water  beyond  the  natural  shoreline.  (This 
apparently  can  help  the  wave  absorption  properties  of  the  beach,  see 
Mahony  a  Pritchard  1980.)  For  the  experiments  to  be  described  the 
forcing  period  ranged  between  0.75  s  and  1.25  s.  The  theory  of 
Mahony  a  Pritchard  (1980)  for  wave  absorption  by  beaches  should  give 
a  fairly  good  estimate  of  the  reflection  properties  of  these  beaches. 
This  theory  indicates  that  a  normally-incident  wave  of  period  0.75  s 
should  have  about  98%  of  its  enery  dissipated  an  the  beach  and  a 
wave  of  period  1.25  s  should  have  approximately  94%  of  its  energy 
absorbed  by  the  beach. 


The  wavenaker 

The  theory  of  Longue t-Higgins  (1967)  indicates  that  the  response 
of  a  given  mode  will  be  abnormally  large  over  only  a  very  small 
bandwidth  of  frequencies.  For  example,  in  the  experiments  to  be  described 
the  period  of  the  wavemotions  was  of  the  order  of  1  s  and  the  'bandwidth' 
of  one  of  the  abnormally  large  responses,  at  about  this  frequency,  is 
predicted  to  be  approximately  0.03  s.  Such  narrow  bandvidths  place 
severe  constraints  on  the  frequency  control  of  the  wavemaker,  requiring 
variations  in  the  wave  period  of,  say,  less  than  3x10  3  s,  a  constraint 
that  has  to  be  maintained  for  a  time  span  of  about  an  hour.  Frequency 
stability  of  this  kind  is  not  usually  possible  with  most  electric  motors 
and  we  therefore  decided  to  base  the  frequency  source  for  the  wavemaker 
on  a  crystal  oscillator  which  was  stable  to  within  5  parts  in  10 7 . 


Over  the  course  of  the  experiment  two  wavenaker  drives  were 


employed.  In  the  early  stages  a  waveeaJcer  of  the  kind  described  in 
Barnard  at  al.  (1977)  was  used,  but  eventually  we  changed  to  a  different, 
■ore  convenient  device.  Pulses  from  the  crystal  oscillator  were 
used  to  control  a  (400  -  pole)  stepping  motor,  the  shaft  of  which  was 
coupled  through  a  crank  to  the  paddle.  The  latter  device  gave  good 
direct  control  over  the  wavemaker  amplitude  as  well  as  providing  the 
frequency  stability  needed  for  the  experiment. 

Measuring  technique 

Measurements  of  the  wave  amplitudes  were  made  using  proximity 
transducers  placed  near  the  surface  of  the  water.  The  principle  of 
the  instrument  is  that  the  transducer  forms  one  plate  of  a  capacitor, 
the  water  surface  being  the  other  plate.  Then,  by  determining  the 
capacitance  it  is  possible  to  infer  the  distance  of  the  water  surface 
from  the  electrode.  The  electronics  associated  with  the  instrument 
(Wayne  Kerr  B731B  Proximity  Meter)  provided  an  output  voltage  that 
varied  linearly  with  the  distance  between  the  two  electrodes,  for 
separations  up  to  about  half  the  diameter  of  the  transducer.  The 
frequency  response  of  the  system  ranged  from  d.c.  to  approximately  1  kHz 

The  output  voltage  from  the  proximity  gauge  was  relayed  to  an 
ultraviolet  chart  recorder,  to  provide  a  permanent  record  of  the 
displacement  of  the  water  surface  at  the  position  of  the  transducer. 

For  a  sinusoidal  response,  which  was  usually  the  case  in  these 
experiments,  the  total  displacement  was  taken  to  be  a  measure  of  twice 
the  wave  amplitude.  The  wave  amplitudes  arising  in  this  study  were  of 
the  order  of  0.05  mm  in  the  main  body  of  the  tank  and  0.25  mm  over  the 
sea  mount.  The  response  over  the  sill  was  found  to  increase  in  direct 


proportion  to  the  amplitude  of  the  wavemaker,  even  at  much  larger 
amplitude  levels  than  those  used  here.  The  capacitance  probes  used 
to  measure  these  waves  had  a  diameter  of  11  mm. 

Measurements  of  the  wavefield  were  made  by  suspending  the  capacitance 
transducers  from  a  moveable  bridge  that  spanned  the  wave  tank.  For  the 
systematic  survey  over  the  island,  seven  transducers  were  mounted  on 
an  arm  that  could  be  rotated  about  a  vertical  axis  aligned  with  the 
centre  of  the  island.  The  transducers  were  located  at  radii  of 
3.6,  10.8,  18.0,  25.2,  32.4,  39.6  and  46.8  cm,  and  measurements  were 
made  over  one  half  of  the  island  at  19  angular  locations,  each 
separated  by  10°  .  (Preliminary  tests  indicated  that  the  wavefield 
over  the  island  was  nearly  symmetric  about  the  diameter  normal  to 
the  wavemaker.) 

.  The  (angular)  displacement  of  the  wavemaker  was  also  determined 
through  the  use  of  a  capacitance  transducer. 

3.  EXPERIMENTAL  RESULTS 

3.1  Notation 

Let  the  origin  for  a  cylindrical-polar  coordinate  system  in  the 
horizontal  plane  be  at  the  centre  of  the  (submerged)  island  and  at 
the  undisturbed  position  of  the  free  surface.  Denote  a  point  x  -  (r,0) 
in  this  plane  by  its  radial  (r)  and  angular  (6)  coordinates,  where 
8  “  O  is  taken  along  the  line  normal  to  the  wavemaker  and  in  the 
direction  from  the  centre  of  the  island  to  the  wavemaker.  Here  r  (and  x) 
are  assumed  to  be  dimensionless  quantities,  being  scaled  in  terms  of 
the  island  radius  a.  Let  the  undisturbed  depth  of  water  over  the 
island  be  d,  and  that  in  the  'ocean'  outside  the  island  be  D.  Define 


the  period  of  the  wavemaker  to  be  T  and  write  0  ■  2w/T.  Let  the 
vertical  displacement  of  the  water  surface  from  its  equilibrium 
position  at  the  point  x  be  u(x,t)  ,  and  define  a  dimensionless  wave 
amplitude  u (xf  t)  -  u(x,t)  /0  ,  where  U  is  a  measure  of  the  magnitude 
of  the  wave train  incident  on  the  island. 

We  suppose  that  u  is  of  the  form 

u(x,t)  -  U(x)exp{i(0(t)  ♦(X))>  , 
where  t  is  time  and  ^  li  •  phase. 

Write  0^  -  U(ri#9j),  where  z±  -  *  j  i  *  i  -  0,1,. ..,6} 

and  ■  {jv/18  t  j  ■  0,1,. ..,18),  and  define  functionals  and  8 
of  U  as  follows) 

0^(0)  -  maxCu^  *  i  -  O,...,  6  i  j  -  0,...,18>  , 

6  18 

and  B(0)  ■  £  £  w;  u*-  r  Ar  A8  , 

i-o  13  X3  1 

where  w^  represents  a  'weight*,  chosen  according  to  some  quadrature 
rule,  associated  with  the  point  (r^,8^)  and  Ar  ■  y  and  A8  ■  w/18. 
Here  the  w^  were  chosen  thus:  for  each  i,  w^^  -  w^^^  -  1, 
wi^2k-1j  -  2  for  k  -  1,....,9  and  w^^  "  4  for  k  (Thaw 

weights  correspond  to  using  Simpson's  rule  as  the  quadrature  formula 
to  approximate  the  8-integration  and  to  using  the  mid-point  rule  to 
approximate  the  integration  over  r.)  The  quantity  Um  is  therefore 
an  approximation  to  the  largest  wave  amplitude  above  the  island 
(at  the  frequency  0)  and  E  is  a  measure  of  the  (square  of  Che) 
energy  in  the  wave field  above  the  island. 

3.2  The  wavefield  over  the  sill 

All  the  experimental  results  described  here  were  made  with 
D  »  15.4  cm  and  d  ■  1.75  cm.  Measurements  of  U.  .  were  obtained  at 


a  number  o£  discrete  values  of  T,  beginning  at  T  -  0.7500  s  and 
incrementing  T  by  0.010  i  to  a  value  of  1.2000  s.  Measurements  of 
U  only  were  then  obtained  with  T  increased  by  the  same  increment 
up  to  a  value  of  1.2500  s.  The  size  of  these  increments  was  chosen 
to  respect  the  bandwidths  estimated  from  ttie  theory  of  Longue t-Higglns 
(1967) . 

(i)  The  amplitude  distribution 


Through  an  interpolation  of  the  measurements  we  have  constructed 
level  sets  for  0  *  0.5  and  0  ■  C  for  C  ■  1,2,3,....  Some  typical 

t 

examples  of  the  contours  are  given  in  figure  2.  All  these  maps  have 
the  common  feature  that  there  is  a  zone  of  relatively  large  amplitudes 
near  the  line  0*0  and  that  the  maximum  wave  amplitudes  were  found 


either  an,  or  very  close  to,  the  line  6  ■  100°.  In  fact,  the  largest 

observed  value  of  0.  .  at  a  given  period  was  either  U or  0,A.  at 

ij  (3)  (18)  (2)  (18) 


all  but  seven  of  the  frequencies  used,  and  in  those  cases  the  a 
response  was  found  at  one  of  the  locations  adjacent  to  these  two.  The 
actual  locations  of  max{0^}  are  given  below  in  table  2.  Further 
discussion  of  the  vavefield  depicted  in  figure  2(e)  is  given  below. 


(ii)  The  functionals  U  and  E 

•  B 

A  graph  depicting  the  measurements  of  0  (T)  and  E(T)  is  given 

n 

in  figure  3.  Also  shown  in  this  figure  are  the  theoretical  values  for 
E(T)  given  by  the  theory  of  Renardy  (1981).  The  reference  wave 
amplitude  used  to  scale  the  experimental  results  in  this  figure  was 
directly  related  to  the  amplitude,  s,  of  the  wavemaker.  Thus,  0  was 


Figure  2a  (T  =  0.8300s) 


Level  sets  for  the  amplitude  field  u(x)  =  u(x)/U. 

(a)  T  =  0.8300s,  U  =  0.233  mm;  (b)  T  =  0.9500s,  U  =  0.226  mm 

(c)  T  =  1.0200s,  U  =  0.284  mm;  (d)  T  =  1.0800s,  U  =  0.139  mn 

(e)  T  =  1.1800s,  U  =  0.208  mm.  The  grid  of  dots  indicate  the 

positions  at  which  measurements  were  made.  The  values  of  u 

are  indicated  near  each  contour.  The  code  is: 


Figure  2c  (T  =  1.0200s) 


Figure  2d  (T  ■  1.0800s) 


Figure  2e  (T  =  1.1800s) 


Renardy  (1981);  the  theoretical  values  for  the  local  maxima  of  U  are  17.5  (near 
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TABLE  2.  Siaairy  of  the  Mia  results.  Ths  quoted  value*  of  UB  end  E  are  those  used  for  figure  3  and  are 

oeasured  relative  to  the  wave  'amplitude*  V  (c  f.  $3.2(11));  W  (c  f.  S3. 4)  is  ths  wave  'amplitude* 
used  to  scale  the  data  for  figure  7#  (I,J)  is  the  location  of  the  sash  point  at  which  aax{u^}  waa 
found.  (Noteei  (1)  0  and  H  are  amplitudes  per  unit  stroke  of  the  waveMkari  (il)  the  eapiricai 
values  for  figurs  7  are  respectively  U/W  and  (U/W) 1  tie* a  the  stated  values  for  U  and  B.) 


taken  to  be  the  (theoretical)  amplitude  of  the  wave fie Id  generated  In 
an  'ocean'  of  uniform  depth  0  by  a  two-dimensional  wavemaker  oscillating 
at  an  amplitude  S,  as  given  by  the  classical  (linear)  theory  of  Havelock 
(1929)  (or  see,  for  example,  Ursell  et  al.  1960).  the  actual  values 
of  U  that  were  used  are  given  below  in  table  2  where  the  main  results 
are  summarized. 

It  is  seen  in  figure  3  that  there  are  a  number  of  similarities 
between  the  empirically  determined  values  of  Om  and  B  with,  for  example, 
their  local  maxima  and  minima  occurring  at  approximately  the  same  values 
of  T.  the  largest  observed  values  of  and  B  were  recorded  at 
T  -  1.1700  and  were  respectively  8.87  and  14.4.  this  means  that  at 
the  most  dangerous  location  above  the  sill  the  waves  were  approximately 
nine  times  as  large  as  would  have  been  expected  with  no  sill  present. 

At  other  periods  the  maximum  amplification  was  not  so  great,  fr*iri«»g 

•  • 

values  of  4  or  S  with  occasional  peaks  and  troughs  outside  this  range. 
The  'base*  value  for  E  was  also  around  4,  (note  that  in  the  absence  of 
the  sill  E  would  take  a  value  of  approximately  ^x) ,  which  is  not 
greatly  different  from  the  'base*  levels  predicted  by  the  theory  of 
Renardy.  On  the  other  hand,  there  is  no  evidence  in  this  set  of  data 
for  the  dramatic  resonances  expected  theoretically. 

The  graphs  shown  in  figure  3  illustrate  the  importance  of  some 
of  the  experimental  considerations  described  above.  For  example, 
the  bandwidth  of  the  'resonance'  near  T  -  0.833s  is  only  about  5o% 
greater  than  the  increment  in  T  used  in  the  experiment  and  the  need  to 
maintain  a  fine  control  over  the  wavemaker  period  for  a  relatively  large 
time  span  is  evident.  An  explanation  of  the  vertical  arrows  labelled 
2,  2.5,  3,  3.5  is  given  below. 
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iii)  The  wavefield  at  T  ■  1.1800a 

Because  the  most  significant  responses  in  this  experiment  would 
appear  to  have  occurred  at  periods  T  between  1.16  and  1.19sf  the 
wavefield  at  T  »  1.1800s  has  been  examined  in  some  detail.  The  data 
from  which  the  contour  map  of  figure  2  (e)  was  constructed  is  depicted 
in  figure  4  where  sections  of  the  field  0(x)  for  constant  r  and 
constant  6  are  given.  TO  help  interpret  these  data  consider  the 
representation  of  the  wavefield  above  the  sill  used  by  Longue t-Higgins 
(1967),  namely 

U(r,0»  t)  -  ZJ^'J'.Otr)  .1<l>e40t)  f 

idiere  the  A^'s  are  complex  constants,  Jw  is  the  Bessel  function  of 
tbs  first  kind  of  order  m  and  y  is  the  real  solution  of 
y  tanh  (dy/a)  ■  aOa/gj  g  is  the  gravity  constant.  It  is  of  interest 
to  find  the  values  for  1,  In  such  a  representation  that  would  yield  ths 
best  fit  with  ths  experimental  data,  through  say  a  Fourier  analysis  of 
ths  data.  But  unfortunately  no  phase  measurements  were  made  in  the 
experiment,  and  a  direct  analysis  was  not  possible.  Nevertheless  an 
attempt  has  been  made  to  'fit'  the  amplitude  data.  Values  for  ths 
(amplitudes)  of  the  coefficients  1,  given  in  table  1  correspond  to 
coefficients  that  give  a  reasonably  good  description  of  the  measurements. 
(The  arguments  leading  to  this  choice  are  rather  complex,  but  they  are 
fully  detailed  in  the  thesis  of  Yamamuro-Renardy  1980).  Also  given  in 
table  1  are  the  expected  values  from  the  Longue t-Higgins  (1967)  theory 
for  a  decomposition  that  realises  an  incident  plane  wave  from  the 
direction  6-0. 


-19- 


Experiment 

1.6 

15.5 

small 

3.1 

small 

Theory 

1.6 

8.0 

3.0 

3.0 

1.6 

Table  1.  Estimates  of  the  magnitudes  of  the  first  four 
FOurier-Bessel  coefficients  Am  that  give  the 
'best*  description  of  the  data  of  figure  4.  The 
theoretical  values  are  those  for  the  theory  of 
Longuet-Higglns  (1967)  at  the  same  frequency. 

This  comparison  suggests  an  almost  complete  absence  of  the 
m  »  2  and  m  ■  4  modes  and  the  enhancement  of  the  mode  m  ■  1.  Attempts 
to  resolve  these  differences  through  nonlinear  effects  mere  unsuccessful , 
though  the  inclusion  of  a  reflected  wave  from  the  beach  behind  the  sill 
(in  the  direction  0  ■  I)  did  improve  the  agreement  a  little.  Thus  it 
would  appear  that  the  resolution  of  these  differences,  especially  the 
large  response  of  the  mode  m  ■  1,  is  unlikely  to  be  achieved  through 
an  analysis  based  on  the  decomposition  of  a  plane  incident  wave. 

3.3  The  wave fie Id  in  the  'ocean* 

Some  measurements  were  made  of  the  wave  field  in  the  region  between 

the  sill  and  the  wavemaker  to  ascertain  the  level  of  any  st»nding-wave 

motion  in  that  zone.  'One  such  measurement  was  made  along  the  line  AB 

indicated  in  figure  1,  wave  amplitudes  being  recorded  at  seven  different 

locations.  Some  examples  of  these  results  are  given  in  figure  5,  where 

the  abscissa  represents  the  distance  from  the  wavemaker  and  0  is  the 

o 

wave  amplitude  relative  to  the  amplitude  at  the  waterline  of  the 
(sinusoidal)  displacement  of  the  wavemaker.  Also  shown  in  the  figure  is  the 
theoretical  wavelength  X  of  the  motion  at  the  given  period.  The  curves 
connecting  the  data  points  were  drawn  only  as  a  possible  guide  to  help 

correlate  the  data,  but  noting  that  with  a  two-dimensional  reflected  wave field. 


D  should  undulate  at  a  wavelength  X/2  (see,  for  example,  Mahony  a 
o 

Pritchard  1980).  (Since  at  the  smaller  values  of  T  the  measurements 

of  0  were  each  separated  by  about  half  a  wavelength  of  the  forced 
o 

motion,  the  dashed  curves  have  been  drawn  just  to  link  the  data  points 
without  respecting  the  anticipated  (X/2)  -  undulations  of  UQ.) 

The  results  shown  in  figure  5  Indicate  the  presence  of  an 
appreciable  standing-wave  component  of  the  wavefield,  arising  as  a  result 
of  reflections  from  the  submerged  cylinder  and  the  beaches.  Since  this 
standing-wave  field  could  easily  give  rise  to  anomalous  results  should  there 

be  any  tuning  effects  (e.g.  between  the  wavemaker  and  the  sill)  we  have 

» 

indicated  by  the  vertical  arrows  in  figure  3  some  of  the  periods  at 
which  such  effects  could  have  been  prominent}  shown  near  each  arrow 
is  the  number  of  wavelengths  between  the  wavemaker  and  the  'front*  edge 
of  the  8 ill.  The  indicated  periods  are  seen  to  coincide  with  local 
maxima  of  t!m  and  E,  except  for  the  one  near  T  =  1.22s  which  is  near  a 
local  minimus. 

A  few  measurements  were  also  made  along  a  line  parallel  to  the 
wavemaker  at  a  distance  56  cm  (=  1.12a)  from  the  wavemaker  to  examine 
the  degree  of  nonuni f ormi ty  in  this  direction.  The  results,  which 
are  given  in  figure  6,  confirm  the  complex,  three-dimensional  nature  of 
the  wave  scattering  from  the  sill  and  the  beaches  already  suggested  by 
the  results  of  figure  5.  Thus,  in  view  of  the  apparently  complicated 
structure  of  the  wavefield  in  the  'ocean'  beyond  the  sill,  it  was  decided 
to  make  an  alternative  assessment  of  the  results  depicted  in  figure  3. 

3.4  A  rescaling  of  the  data 

A  possible  explanation  of  the  result  that  none  of  the  theoretical 
resonances  were  apparent  in  the  results  of  figure  3  is  that,  because  of 


-23- 


Wave  amplitudes  measured  along  a  line  parallel  to  the  wavemaker  at  a  distance  1.12a 
from  the  wavemaker.  y  is  the  distance  from  the  centreline  of  the  tank;  the  amplitude 
scale  is  arbitrary;  A  and  B  represent  the  positions  of  the  left  and  right-hand 
sides  of  the  wavemaker.  - :  T  ®  0.8015s;  — — — — :  T=  1. 0015s ;  • - :  T=  1.2015s 


the  appreciable  standing-wave  component  of  the  field  in  the  'ocean' , 
the  'estimate'  for  the  effective  incident  wave  used  there  was  grossly 
inaccurate.  We  have,  therefore,  tried  to  assess  the  importance  of 
this  feature  by  using  data  of  the  kind  shown  in  figure  5  to  provide  an 
empirical  estimate  W,  say,  of  the  incident  wave  field  from  the  'mean' 
wave  amplitude  in  the  "zone  between  the  wavemaker  and  the  sill  (c  f. 
Mahony  ft  Pritchard  1980  for  a  description  of  this  kind  of  partitioning). 
For  example,  at  period  T  -  1.22s  for  the  data  shown  in  figure  5(f), 

W  was  chosen  to  be  0.174.  All  the  values  for  W  thus  chosen  are  given 
in  table  2,  together  with  the  values  of  U  used  to  construct  figure  3. 

Die  result  of  using  the  estimate  W  for  the  incident  wave  levels  is  shown 
in  figure  7,  where  and  E  have  been  replotted.  Here  the  four  largest 
maxima,  in  order  of  decreasing  magnitude,  occurred  for  E  at 
T  -  1.19,  0.96,  0.80  and  1.08s  and  for  at  T  -  1.18,  0.96,  1.08  and 
0. 80s.  Note  that  the  overall  maximum  for  figure  3  was  near  T  -  1.17s, 
so  that  the  rescaling  has  not  seen  the  elimination  of  the  large  response 
near  this  period.  But  while  some  more  prominant  local  maxima  are 
evident  in  figure  7,  there  is  still  no  indication  of  the  very  large 
responses  expected  theoretically. 

4.  ASSESSMENT  ft  DISCUSSION 

The  results  obtained  from  these  experiments  do  not  show  any 
indication  of  the  resonances  expected  theoretically  when  a  train  of 
plane  waves  passes  over  a  submerged  cylindrical  island  situated  in  an 
otherwise  uniform  oceah.  Some  of  the  possible  reasons  for  this 
disparity  have  been  examined  in  S3  where,  in  particular,  attempts 
have  been  made  to  assess  the  laboratory  experiment  in  terms  of  the 
hypothetical  problem  studied  theoretically. 
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The  elimination  of  unwanted  motions  in  water  wave  experiments 
is  usually  very  difficult  (e.g.  see  Bona,  Pritchard  ft  Scott,  1981)  and 
the  present  experiment  has  proved  to  be  no  exception.  Thus,  a 
particular  difficulty  arising  in  this  study  has  been  that  of  defining 
a  reference  wave  amplitude  against  which  the  response  over  the  sill 
could  be  scaled.  Two  reference  amplitudes  were  considered  in  13  : 
one  is  related  directly  to  the  amplitude  of  the  wavemaker  used  to 
force  the  motions  (through  classical  wavemaker  theory)  and  the  other 
was  determined  from  wave  measurements  in  the  tank.  Bach  of  these 
procedures  is  inadequate  in  one  way  or  another  (with  regard  to  their 
relation  to  the  theoretical  situation) .  Nevertheless ,  at  periods 
remote  from  the  theoretical  'resonances',  the  measured  values  of  tile 
quantities  UM  and  E  (c.f.  figures  3,  7)  were  in  reasonably  good 
agreement  with  the  theoretical  predictions  of  the  model  calculation 
in  the  sense  that  the  discrepancies  between  the  theory  and  the 
experiment  could  be  accounted  for  by  a  change  of  the  reference-wave 
amplitude  of  not  more  than  about  30%.  Near  the  resonant  periods, 
however,  the  estimates  for  the  reference-wave  amplitude  would  have  to 
be  changed  by  a  factor  of  over  3  to  account  for  the  observed  discrepancies, 
and  such  large  factors  should  have  been  apparent  in  the  experimental 
results.  An  important  aspect  of  the  discussion  of  S3,  and  bearing  on 
the  above  issues,  was  that  the  response  over  the  sill  was  not  entirely 
compatible  with  the  decomposition  of  an  incident  plane  wave  from 
Infinity  and  this  would  appear  to  present  the  main  source  of  difficulty 
in  interpreting  the  experimental  situation  in  terms  of  the  hypothetical 
model. 


A  salient  feature  of  the  'resonance'  phenomenon  wider  consideration 
is  that  leakage  rates  of  certain  eigenvalues  become  small  near  particular 
frequencies,  enabling  the  build-up  of  the  corresponding  wavemode  to 
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unusually  largu  amplitudes.  But,  in  consequence  of  the  saall  leakage 
rate#  the  coupling  between  the  external,  incident  wave  fie  Id  and  that 
over  the  island  is  very  weak#  so#  for  significant  amounts  of  energy  to 
be  transferred  to  the  field  above  the  island#  the  exterior  conditions 
must  be-  held  steady  for  a  considerable  time. 

The  design  of  the  wavemaker  drive#  together  with  the  very  close 
control  maintained  over  the  frequency  input  to  it#  ensured  that 
phase  drifting  was  unimportant  here.  Also#  the  very  small  amplitudes 
(and  surface  slopes)  employed  in  this  study  would  appear  to  have 
precluded  the  possibility  that  nonlinear  effects  could  have  overwhelmed 
the  resonance.  On  the  other  hand#  it  is  more  difficult  to  assess  the 
importance  of  surface  contamination  on  the  motions.  He  do  know  that 
if  the  surface  became  unduly  contaminated  the  wave  fie  Id  was  influenced 
significantly;  but  the  fact  that  the  general  level  of  response  was  in 
broad  agreement  with  the  levels  anticipated  theoretically  (c.f  figures  3,7) 
suggests  that  the  effects  of  surface  contamination  was  not  a  dominating 
feature  of  the  experiment.  According  to  Renardy  (1981)  dissipative 
effects  on  the  rigid  boundaries  should  not  have  been  important. 

Thus#  although  the  experiment  did  not  apparently  realise  a  good 
representation  of  the  hypothetical  wave  fie  Id  needed  for  the  theoretical 
■odel,  it  is  non  the  less  surprising  that  there  were  no  manifestations 
whatsoever  of  the  dramatic  resonances  predicted  theoretically#  and  an 
explanation  of  the  observed  response  would  therefore  be  desirable. 

Because  the  beaches  used  in  the  experiments  should  have  been  fairly 
efficient  absorbers  of  energy#  such  an  explanation  might,  in  the  first 
Pl*ce#  be  sought  from  a  model  with  a  long  wavemaker  generating  waves 
ecu:  an  island  in  a  semi-infinite  ocean.  He  believe  this  experiment 
might  present  a  reasonably  critical  test  for  numerical  schemes  for 


water-wav*  problems  of  this  kind  and  we  have  therefor*  included  (see 
table  2)  a  summary  of  the  sain  experimental  results. 

Finally  we  should  like  to  make  some  observations  concerning  the 
siting  of  structures  on  submerged  ridges  of  this  kind.  While,  in  the 
above  study,  the  theory  indicates  that  at  certain  'resonant'  frequencies 
unusually  large  wave  responses  can  prevail,  both  the  theory  and  the 
experiments  indicate  that,  at  certain  dangerous  locations  on  the  sill  the 
waves  are  magnified  by  a  factor  of  about  five  tines  that  of  the  prevailing 
ocean  swell,  irrespective  of  the  frequency.  Under  moderate  conditions 
at  sea  this  might  not  be  a  serious  problem,  but  with  a  heavy  swell  running 
an  amplification  by  five  times  could  be  catastrophic.  Fortunately  the 
very  dangerous  locations  seem  to  be  fairly  insensitive  to  frequency,  for 
a  given  prevailing  direction)  by  contrast,  there  are  apparently  also 
quiescent  regions,  again  fairly  insensitive  to  frequency,  where  the 
amplication  is  less  than  one. 
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